Intrinsic atomic scale modulations of the superconducting gap of 2H-NbSe 2 
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We present scanning tunneling microscopy and spectroscopy measurements at lOOmK in the su- 
perconducting material 2H-NbSe2 that show well defined features in the superconducting density 
of states changing in a pattern closely following atomic periodicity. Our experiment demonstrates 
that the intrinsic superconducting density of states can show atomic size modulations, which reflect 
the reciprocal space structure of the superconducting gap. In particular we obtain that the super- 
conducting gap of 2H-NbSe2 has six fold modulated components at 0.75 mV and 1.2 mV. Moreover, 
we also find related atomic size modulations inside vortices, demonstrating that the much discussed 
star shape vortex structure produced by localized states inside the vortex cores, has a, hitherto 
undetected, superposed atomic size modulation. The tip substrate interaction in an anisotropic 
superconductor has been calculated, giving position dependent changes related to the observed gap 
anisotropy. 

PACS numbers: 74.25. Jb,74.50+r,74.70Ad 

I. INTRODUCTION 

Today one of the more creative spaces in research related to superconductivity are atomic scale investigations 
by means of nanoprobes on well characterized surfaces of superconductors. Nanometric size inhomogeneities of 
the superconducting local density of states (LDOS) found using scanning tunneling microscopy and spectroscopy 
(STM/S) in high crit ical t emperature superconductors, related to either the pseudogap state of the normal phase, 
or to scattering effect s^* 3 -*, have clearly shed new light on the very nature of the pairing mechanism. However, it is 
generally considered that the superconducting LDOS is homogeneous at nanometric length scales, in ordered, metallic 
and pure samples^ 5 -*. Here we report on the discovery of intrinsic atomic scale modulations of the superconducting 
gap using STM/S measurements at 100 mK. 

We have found these modulations in one of the superconducting materials most commonly used in STM/S studies at 
low temperatures, the layered dichalchogenide superconductor 2H-NbSe 2 . This compound has a critical temperature 
of T c =7.1 K and develops a charge density wave (CDW) state below T CDW =33 K. Historically, it has been one of the 
first superconducting materials where the superconducting LDOS has been studied in depth. This is not surprising, 
as it is easy to prepare flat and clean surfaces. Indeed, 2H-NbSe2 is composed by layers of one sheet of Nb atoms 
between two sheets of Se atoms (Se-Nb-Se). Adjacent layers are bonded with weak Van der Waals forces which leads 
to a pronounced two-dim ensio nal character, making this material very easy to exfoliate. The last layer (surface layer) 
is formed by the Se atoms^*^*. First observations of the vortex lattice with STM/S in a superconductor were made in 
this materiaP^ 1 . It is now believed that 2H-NbSe2 is an important system where effects concerning interplay between 
charge order, superconductivity and a non trivial Fermi surface structure can be handled and studied in detail. 
The deviations from simple s-wave BCS theory of many experimental properties in this compound are notable and 
provide interesting analogies t o supercondu ctors close to so called quantum critical instabilities, and to the problem 
of multiband superconductivit y 11 * 12 * 13 * 14 * 15 ! The Fermi surface consists of several nearly two dimensional concentric 
cylinders, de rived from the Nb 4d orbitals and a single small three dimensional pancake like sheet, derived from the 
Se 4p orbitalsPE2. The local superconduc ting de nsity of states shows quasiparticle peaks far from the high anomalies 
expected from simple isotropic BCS theory^Mi8*. This evidences the presence of a wide spread distribution of values of 
the superconducting gap, with gap sizes ra nging between 0.6 meV and 1.4 meV, i.e. several hundreds ^eV around the 
value expected from BCS theory (1.1 meV^Esl Further experiments highlight a close relationship between the band 
structure, the distribution of values of the superconducting gap, and the CDW state. The wide spread distribution of 
values of the superconducting gap have been proposed to reflect multi-band superconducting properties, with different 
superconducting gaps in th e different Fermi surface sheets of this compound, in analogy to the intriguing two band 
superconductor MgBjp^H. In particular, thermal conductivity and specific heat studies have shown the presence 
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of at least two well differentiated values of the superconducting gap over the Fermi surface^ESl. Angular resolved 
photoemission experiments (ARPES) have been made only relatively close to T c (at 5.3 KjP^I, to find that Nb band 
gaps are around 1 meV. However, the Se band showed no gap opening, although London penetration depth studies 
give strong evidence for a relatively large gap in this band at low temperatures^. Moreover, de Haas van Alphcn 
experiments in the mixed state and at very low temperatures detected a small superconducting gap of 0.6 meV, which 
was associated to the Se sheetP^. Tunneling spectroscopy experiments demonstrated that the smaller size gap feature 
closes indeed rather fast (around 5.5 K) when increasing temperature^. On the other hand, most recent ARPES 
experiments^! made with beautiful detail (although also at 5.3 K, which corresponds to state of the art experimental 
possibilities), have shown that the maxima of the superconducting and CDW gaps appear at the same position in 
reciprocal space, suggesting a cooperative interaction between charge fluctuations and superconductivity. This appears 
to go against the view obtained from pressure experiments, where it is clear that application of pressure decreases 
Tcdw, while T c increases^. However, T c continues increasing after the destruction of the CDW state, and the upper 
critical field shows also a peculiar and complex pressure dependence^, demonstrating that the pressure induced 
changes in the Fermi surface are non trivial and could produce in some unknown way a simultaneous destruction of 
CDW and an increase in T c . On the other hand, the peculiar form of the internal electronic structure of vortices 
has al so b een put in connection with the ongoing debate about the anisotropic and/or multiband properties of 2H- 
NbSegpnS Experiments show that the LDOS in and around vortices has a characteristic star shape with a complex 
energy dependence, that has been explained in a great deal by taking into account the symmetry of the vortex lattice, 
the formation of quasiparticle bound states inside vortex cores, as well as a six-fold modulated superconducting gap^l 
However, the connection between this six-fold modulation and the zero field gap distribution, as revealed from the 
experiments, has not been made until now. 

Within this debate, it is of importance to provide new precise information about the anisotropy and values of the 
superconducting gap at the lowest possible temperatures, i.e. close to the superconducting "ground state". Indeed, 
in a multiband superconductor the temperature dependence of the gaps found are expected to be non-trivial. Strong 
evidences for this has been obtained from the mentioned tunneling spectroscopy experiments, where changes in the 
distribution of values of the superconducting gaps in 2H-NbSe2 have been followed as a function of temperature^!. 
Moreover, it is equally important to relate the mixed phase observations with the zero field results. The atomic scale 
modulations of the superconducting gap reported here (found below 0.015T C ) show directly that the superconducting 
density of states has sixfold modulations at 0.75 mV and 1.2 mV, and therefore demonstrate that the associated 
gaps are six-fold modulated in reciprocal space. Measurements made under magnetic fields close to the vortex core 
demonstrate that these six fold atomic size modulations of the superconducting gap are superposed to the previously 
detected features of the LDOS inside the vortices. 



II. EXPERIMENTAL 

We use a very low temperature STM/S system installed in a dilution refrigera tor. Th e set-up allows for in-situ 
clean preparation of tip and sample, following procedures described somewhere els d 28 ! 29 * 3 -^. To ensure high resolution 
in voltage, all wires are filtered or attenuated. The experimental setup has been tested thoroughly before through the 
measurement of several low critical temperature superconductors such as Al (T c = 1.12 K), where we obtain clean BCS 
s-wave tunneling conduc tance curves with a negligible conductance at zero bias. The resolution in energy of our set-up 
lies below 15/neV ^ 29 l 31 [ The atomic resolution topographic (STM) images are taken at constant tunneling current. 
The Se hexagonal lattice a nd the charge density wave are observed very clearly in 2H-NbSe2 as in previous STM work 
on this compound (see e.g! 7 * 10 * 31 ^). Full tunneling conductance vs. bias voltage a{V) curve ranging from voltages well 
above the superconducting gap are taken at each point of the topographic images, to construct STS images. The 
place where the cr(V) curves are taken is uniquely identified with the atomic positions through comparison with the 
topography. The whole set-up is isolated from acoustic noise as well as vibrations from pumps and the building to 
allow for long term measurements. As a result, the conductance can be measured at each position with an accuracy of 
less than several per cent. In order to show the effects discussed below, we have to make a very close scan at atomic 
level. About a hundred tunneling conductance curves are taken for each group of three surface atoms, and full images 
cover regions including up to fifty x fifty atomic rows. 

III. RESULTS AND DISCUSSION 

All a(V) curves are very similar to each other, and to the curves obtained in previous measurementJ^^l. The 
tunneling conductance smoothly increases at about 0.5 mV, has a small shoulder at about 0.75 mV and a rounded 
peak around 1.2 mV. However, here we are able to find slight, but relevant, differences between the curves taken at 
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FIG. 1: (Color online) In the lower right inset of a we show typical topographic images obtained in 2H-NbSe2, where the atomic 
Se lattice is clearly viewed, as well as the charge density wave, which produces an additional modulation in topography each 
three Se atoms. In the figure we show a(V) curves obtained on top of a Se atom and in between. STS images showing small 
groups of atoms constructed from these data at Vi = ± 0.75 mV are shown in the top panels. In b we show da / dV (V) at the 
same positions. The panels in the figure show STS images constructed from da/dV at Vi. Full range in contrast in the STS 
panels from black to white represents, respectively a difference in a of 0.3/i S (a) and in da/dV of 2mS/V (b). Circles in the 
images give the Se atomic positions as determined through topography. 



the top of the Se atoms, and those taken in between. In Fig.[T]a we show tunneling conductance curves taken at these 
two positions. The shoulder at Vi = ± 0.75 mV is better defined on a(V) taken over a Se atom. Correspondingly, 
STS images constructed from these curves at Vi (top panels of Fig.[T|a) clearly show a modulation of the conductance 
as a function of the position which reproduces the lattice of Se atoms. Note that this modulation corresponds to small 
changes in the tunneling conductance of the order of 10%. The derivative of <r(V) (Fig. [l]b) has a peak around Vi 
and another one at V2 = ± 1.2 mV. STS images obtained from the derivative of a(V) at these voltages also clearly 
show the surface atomic lattice, as presented in the panels of Fig.[T]b, where we show the result at Vi. 
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FIG. 2: (Color online) In the inset we show the Fourier transform of an STS image constructed from a(V = 0.75mV) in 
2H-NbSe2 at 100 mK. Red circles represent the position of peaks in the Fourier transform of the corresponding topography 
image. The Fourier transform of STS images at other voltages is either flat, close to the Fermi level and at high voltages, or 
has the same form as the one shown. This is seen in the figure, where we represent the height of the six fold peaks of the 
Fourier transform shown in the inset as a function of the bias voltage (line is a guide to the eye). It is important to remark 
that no other peaks appear in the Fourier transform when changing the voltage, and that the features discussed in this paper 
fully disappear at voltages well above the superconducting gap. 

Actually, the full bias voltage dependence of the size of the conductance modulation is best discussed by analyzing 
the Fourier transform of the STS images. In the inset of Fig. [2] we show the result of the Fourier transform of STS 
images obtained from the conductance at Vi. We observe a clear sixfold structure with Fourier peaks whose position 
coincides with the six fold peaks due to the surface Se lattice obtained in the Fourier transform of the associated STM 
topographic image. The height of these peaks strongly changes with bias voltage (Fig. Ejl. The six fold modulation 
appears around 0.6 mV, becomes best denned at Vj,, with again a small feature at V2, and disappears for higher 
voltages. 

Note that, contrasting the results found in other superconducting systems, related to impurities and imperfections 
at or near the surface^, the position of the Fourier peaks we find here does not change with voltage, i.e., there is 
no dispersion with energy. Indeed, periodic structures that have a significant dispersion in energy can be related to 
surface wave function scattering effects, which are much discussed in simple metals and in High T c superconductors, 
as the y ma ke it possible to obtain the electron dispersion relation by analyzing the Fourier transform of the STS 
images^m. This is not the subject of present results, where we focus on measurements on neat surfaces, with a well 
ordered Se lattice. Here, no additional periodic structures, other than the atomic lattice, are resolved in the STS 
images at the voltage range studied here. These periodic structures are found only well within the quasiparticle peaks, 
and not at voltages above the superconducting gap, and show no changes of its reciprocal space position as a function 
of the bias voltage. 

In the mixed state we also find an atomic scale modulation of the superconducting features. It is well known that, 
close to the core of isolated vortices, the superconducting LDOS of 2H-NbSe2 s hows a complex position dependence at 
length scales comparable to the superconducting coherence length £ sa 10 nirpUSI. The whole voltage dependence of 
the LDOS of the vortex lattice of 2H-NbSe2, as found here in full agreement with Ref. 10 , is given in EPAPS Document 
No. []. Time scale reproduces the bias voltage of the STS experiment. The video begins at Vsias = - 2.5 mV and 
ends at V Bias + 2.5 mV, and each sequence corresponds to a change in bias voltage of 0.1 mV. The star shape shown 
in the inset of Fig|3]a is found close to the Fermi level, which appears at the middle of the video. Indeed, the way 
the gap closes inside a vortex core is not trivial, and one of the most characteristic features is the six fold star shape 
structure which is observed on STS images obtaine d from the zero bias conductance. In the inset of Fig. [3] a we show 
the star shape of a single vortex at 0.03 T (see also 1 27 * 32 *). 

Close to the center of the vortex, at the small square represented in the inset of Fig. [3] a, STS images made at 
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FIG. 3: In the inset of a (lower left corner) we show the STS image made at zero bias of a single vortex at 0.03 T. Color scale 
represents the value of the LDOS at the Fermi level, following previous representations of a vortex as measured with 
(the whole color change corresponds to differences in the conductance from zero to about 2 /jS). A zoom up of this image close 
to the vortex core, at the small red square of the inset, gives the STS images a-e, made at the bias voltages represented in 
each panel. Color scale in these panels, from green to blue, represents differences in the conductance of 0.2 /j,S, i.e. changes 
not higher than 10% of the high bias conductance. The bias dependence of the height of the Fourier transform peaks of the 
images is shown in f. Points which have been obtained from the Fourier transform of images a-e are highlighted as greater 
circles. The inset shows the Fourier transform of the STS image at 0.2 mV, with red circles signalling the position of the peaks 
from the Fourier transform of the topography. Again, the effect of the bias voltage is to change the height of the peaks shown 
in the inset, becoming featureless at high bias voltages, as shown by the decrease of the Fourier amplitude in the figure. No 
new significantly remarked peaks appear when changing the voltage. 



atomic length scales again show the atomic surface lattice, only at voltages where superconducting features appear 
on the LDOS. In Figs. [3] a-e we show STS images at different bias voltages. Note that the zero bias voltage image 
shows a smooth color change from more blue at the bottom to more green at the top. There is a smooth increase in 
the conductance, superposed to the atomic size modulation, which reproduces the smooth variation in the tunneling 
conductance that leads to the star shape shown in the inset of Fig(3] a. This variation is less pronounced at higher 
bias voltages, although it is always present and superposed to the atomic size modulations discovered here (Figs. [3] 
a-e) . Note also that the star shape structure is aligned with the hexagonal modulation of the superconducting gap 
and with the atomic lattice. The bias voltage dependence of the height of the corresponding Fourier peaks (Fig. [3]f) 
is very different than the one found at zero field. Although in both cases no modulation is observed at voltages well 
above the quasiparticle peaks, the modulation under magnetic fields is observed also close to the Fermi level, as the 
superconducting LDOS is not zero close to the vortex core. A slight super-modulation at the reciprocal wavevectors 
of the charge density wave is also observed in some of these images. The height of the corresponding Fourier peaks 
are much smaller than the ones shown in Fig.|3]f, and its energy dependence is very difficult to extract from present 
data, although they appear to be better defined close to the Fermi level. Higher resolution experiments may be able 
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to resolve these details in future. 

It is important to remark that the atomic surface lattice is viewed in the STS images only at voltages of the 
order of the position of the superconducting features in the LDOS. For voltages well above the highest value of the 
superconducting gap, the STS images are featureless, as seen in the strong decrease in the height of the Fourier peaks, 
shown in Fig.[2]and in Fig.[3]f. This rules out any effects related to tip-substrate coupling phenomena not intrinsically 
due to the local superconducting properties. 

It is well known that the atomic scale LDOS of the sample, as viewed in the STS experiment, is formed by 
contributions from different parts of the Fermi surface, weighted by details of the interaction between tip and sampleP. 
Many approaches are being developed and applied to understand spectroscopic imaging in different kinds of metallic 
or semiconducting surfaces. The superconducting state has received less attention, probably because it is generally 
considered that, in a first approx imation and in pure systems, the superconducting gap is relatively homogeneous over 
distances of the order of ^ J 4 I 5 I 33 | Nevertheless, it is clear that non trivial variations of the superconducting density of 
states in reciprocal space, as found in multiband superconductors or in systems with a wide distribution of values of 
the superconducting gap around the Fermi surface, have to appear as modulated structures in the real space LDOS. 
Therefore, the modulated features observed on the STS images at Vi and V2 show the presence of concomitant 
modulations of the superconducting gap on reciprocal space. The six fold modulations at Vi = 0.75 mV and V2 = 
1.2 mV observed in our experiments are associated to sixfold reciprocal space modulations of the superconducting 
gap of 2H-NbSe 2 . 

The Fermi surface position where this six fold modulated gap structure appears is more difficult to establish 
independently from our experiment at present. In the model described below, we give a first approach to extract more 
detailed reciprocal space information of the superconducting gap from its atomic size modulations. It is however useful 
to first compare our result with previous experiments. For instance, while the size of the distribution of values of the 
superconduct ing gap ha s been long known, since the seminal tunneling experiments of Hess et alJ 9 * 10 ! and subsequent 
measurements^- 19 * 22 * 23 !, the discussion about the reciprocal space dependence has been limited to the values of the 
different gaps in different sheets of the Fermi surface. Eventual anisotropies and its symmetry remained difficult to 
extract directly from the experimental data available. For example, the connection of the gap anisotropy and the 
six fold star shape structure of the vortex cores, found in RefpSl was made through the ad-hoc introduction of a 
modulated superconducting gap in theoretical calculations^. 

On the other hand, very recent ARPES experiments have shown that the gap has maxima at points in reciprocal 
space that are directly connected with the CDW vector 2 ^. A direct comparison to our experiment is difficult, because 
ARPES measurements were made, although with great detail, at temperatures very close to T c (at 5.3 K, while 
T c =7.1 K), and the fits to the spectra obtained needed the addition of a smearing parameter to get the values of 
the superconducting gap. Our data, by contrast, are taken at 0.1 K, and have a flat and zero conductance at the 
Fermi level, giving neat spectra that do not show any sign of smearing effects. The ARPES experiment shows that 
there are six-fold modulations of the superconducting gap along directions that coincide with the atomic lattice and 
CDW reciprocal wavevectors. Some of the corresponding gap magnitudes are very similar to the ones found here (e.g. 
1.2 mV and 1 mV along the T-K symmetry line 2 ^). Indeed, the anisotropy observed in the present experiment is 
very close to the one reported in some of the ARPES results. Although the CDW wave-vectors are practically not 
observed in the Fourier analysis of the atomic size modulations of the superconducting gap shown here (except in the 
experiment shown in Figs. [3] where the corresponding signal is very weak), the directions of the sixfold anisotropy 
observed here is the same as the CDW wave vector direction. 

Additional angular changes of the superconducting gap are also found in ARPES data for significantly lower values 
of the gap, going down to 0.3 mV. This is fully absent in our experiment, as the LDOS, at zero field, is flat and equal 
to zero (within one pe r cent of the high voltage value) below 0.5 mV. However, as mentioned previously, the tunneling 
experiments of RefFSS demonstrate that the smaller sized gaps of 2H-NbSe2 close when approaching T c , so that the 
apparent discrepancy is clearly due to temperature differences. Note that therefore, a significantly sized anisotropy 
in the Se sheets at low temperatures, which may have an additional influence to the results found here, should not be 
fully excluded at present. 

Anyhow, our experiment gives a determination of the values of the six-fold anisotropy in the superconducting gap 
with unprecedented detail, free from temperature and smearing effects, and also additional local information close to 
the vortex cores. Indeed, the atomic size modulations found within the vortex core should come as a consequence, 
within these considerations, of the closing of the superconducting gap, and produce the changes in the height of the 
Fourier peaks towards energy values closer to the Fermi level. These changes, however, are largely non trivial, as seen 
from Fig. [3] f. First, the Fourier peaks shown in Fig.[3]f reveal that the smaller sized gaps tend to close earlier when 
approaching the core, so that the anisotropic gap at V2 remains stable pretty close to the core, and the anisotropic 
gap at Vi shifts to the Fermi level. Indeed, the Fourier peak at V2 remains practically at the same position as in zero 
field. And second, the sixfold modulation of the superconducting gaps is maintained in its full importance well within 
the vortex core, as evidenced by the height of the Fourier peaks, which are of the same order at V2 than close to the 
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Fermi level, e.g. at Vi (Fig. [3]f). So the effect of the magnetic field inside the core is not a simple smearing effect, 
and the much more tricky formation of localized states, combined with gap anisotropies, as introduced in Ref.^-M, must 
be considered as a central issue in the interpretation of vortex lattice images using STS. 



IV. THEORY 



The development of a microscopic model of the dependence of the observation of superconducting properties as a 
function of atomic scale tip substrate interactions, taking precisely into account the whole Fermi surface of 2H-NbSe2 
is out of the scope of the present paper. Nevertheless, we have verified that a simple description of the band structure 
and tip substrate interaction in 2H-NbSe2 leads to intrinsic modulations of the superconducting density of states 
on atomic length scales, related to the Fermi surface dependence of the superconducting properties. Within this 
description, the anisotropies of the order parameter lead to modulations in the density of states measured by a local 
probe on length scales comparable to the lattice spacing. The two ingredients of the model which lead to changes in 
the superconducting state are the degeneracies of the electronic bands at high symmetry points in the Brillouin Zone, 
and a modulation of the superconducting gap, compatible with the lattice and with previous proposals to explain 
some features of 2H-NbSe2 . 

First, we present the model Hamiltonian we use. We approximate the band structure by a nearest neighbor tight 
binding model which includes the d x i_ y i and d xy orbitals. We consider a single hexagonal Nb layer. This model is 
an approximation to a more complex three dimensional model which includes Se p z orbitals and Nb s and the other 
Nb d orbitals. If the interlayer couplings are weak, and the Se and Nb orbitals not included initially in the model 
are sufficiently far from the Fermi surface, they can be included perturbatively. The resulting hamiltonian is given 
by a 2 x 2 matrix, with the same degeneracies and symmetries as one defined solely in terms of the d x i_ y i and d xy 
orbitals. 

The model includes two nearest neighbor hoppings, tdda and tdd-n- We define: 
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where a is the lattice constant. The hamiltonian is: 
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We can expand the hamiltonian around k x , k y = 0: 
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which can also be written in angular coordinates as: 
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where X is the (2 x 2) unit matrix, and a z and a x are Pauli matrices. The eigenvalues of this matrix are 



k 2 a 2 
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and the eigenfunctions are: 
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where A: is a vector at the Fermi surface, 9 T is the angle which denotes the position of k ar the Fermi surface and \ddn) £ 
nd \dda)j: are Bloch waves built from ddn and dda atomic orbit als. The bands are isotropic, and the eigenvectors 
rotate around the T point, with a winding number of 2n. 

A similar expansion can be done around the K and K' points. We define, for the K point: 



7 n r, 

k x a = — + dk x a 

k y a — — h 5k y a (7) 

o 

The hamiltonian becomes: 
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with eigenstates 

and eigenfunctions: 



|%) = - sin f^j I ddir) g + cos — J |ddo-)g (11) 

The winding number in this case is 7r. 

The above analysis justifies the almost circular Fermi surface sheets obtained from diagonalizing the hamiltonian, 
eq.[2], shown in Fig. [4] a. The symmetries of the eigenfunctions are also consistent with the general behavior of a 
Hamiltonian with twofold degeneracy defined in an hexagonal lattice. In particular, the linear (Dirac) dispersion 
found near the K and K' points is analogous to that observed at the corners of the Brillouin Zone in graphene 34 35 . 

The model can be enlarged by including the d z i orbitals of Nb and the p z orbitals of Se, and also interlayer hoppings, 
leading to four bands which disperse in all directions. For reasonable parameters, however, the results do not change 
qualitatively with respect to the simpler two band model described here. 

We now consider the tip-substrate interaction. The current through the tip measures the LDOS. If the tip is not 
placed on top of a Nb site, it will couple more strongly to a certain combination of d x 2_ y 2 and d xy orbitals. We 
consider a tip which interacts with more than one atom. The interaction is assumed to be anisotropic. In Fig. [4] 
b we show schematically this anisotropy through an ellipse centered at the position of the tip. We assume that the 
tunneling current is determined by a single combination of atomic orbitals at the tip, \tip). The current will be 
proportional to the overlap of this combination of orbitals and the atomic orbitals at the substrate. These overlaps 
can be written in terms of a and 7r components. 

The relevant orbitals in the substrate are Bloch waves near the Fermi energy, derived from the atomic orbitals 
as discussed previously. We consider the Fermi surface sheet near the K point, and the Bloch waves \d x 2_ y 2)^ and 



\dxy)ti 33 defined in eq.( 11 1. A STM tip, at some distance of the Nb atoms, will couple to a combination of the d x 2_ y 2 
and d xy orbitals, whose in plane orientation with respect to the lattice axis can be described by 9 . 9 is the angle 
between the long axis describing the anisotropy of the tip orbitals (long axis of the schematic ellipse shown in Fig. [4] 
a), and the line between the tip and the nearest located atom. Then, the overlaps are, for the a contribution: 
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FIG. 4: In a. we show Fermi surface obtained from the model discussed here (a is the lattice spacing). High symmetry points 
are labeled as F, K and K' . In b. we represent schematically the anisotropy of the tip orbitals. The ellipse in the figure 
graphically shows the anisotropy of the tip-substrate interaction. The Nb d x 2_ y 2 and d xy orbitals are represented as in plane 
lobes (green and red respectively). 80 is the angle between the long axis of the anisotropy of the tip and the straightest line 
connecting the center of the tip to the nearest Nb atom, schematically represented here at the center of the 2D real space 
coordinate system. 



where the labels i,j stand for lattice positions, defined with respect to the axis shown in Fig. [5] and the angular 
coordinates , 0y are used. There is a radial modulation of the contribution from each lattice site, included in the 
function f(rij). The function J%(x) is a Bessel function. 9 K is used for the angular coordinate in reciprocal space with 
respect to the K point, as previously. A similar expression is found for the 7r contribution, except that the factors 
cos[2(0 K — 9q)] and sin[2(# K — 9q)] are interchanged. 

The anisotropic contributions in eq.(12l have to be integrated over the Fermi surface. As a further ingredient, 
we need to give the simplest non isotropic modulation of the superconducting gap , compati ble with the hexagonal 
symmetry of the lattice (extended s-wave with a finite non-isotropic component, se d 12 l 26 * 2 ^^j), which can be written 
as Ag + cos(30 ) around the K point. Then, the density of states measured by the tip includes a position 
dependent term: 



de 



K 



cos[2(0 K - O )] cos + sin[2(0 K - 0„)] sin (e£\ 

y/u; 2 - [Af + Af cos(36^)] 2 



(13) 



A similar equation describes the contribution from the K 1 point, with A^ = — Af^. The hexagonal symmetry is 
recovered from adding both terms. 

The same analysis, when applied to the states around T leads to the correction to the Green's function: 
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FIG. 5: (Color online) In the figure we plot the LDOS for different angles when the tip describes a circle around a Nb 
atom (do = 0°, black, 10° in red; and 20° in green). We show the result near the K point. The parameters used for the 
superconducting gap are Ao = 1 meV and Ai = 0.3 meV. 



6G Up ,r(u) = / d9 



|cos[2((9 r - 9 )] cos(0 r ) + sin[2((9 r - 6 )] sin(6» r ) f u 



[A 



3 (60 r )] 2 



(14) 



The expression for 5Gti P ,K in eq.(13l leads to a significant dependence on 9q, as the numerator when expanded in 
angular harmonics, includes components proportional to W K , as in the denominator. This is not the case for the 
contribution from the Fermi surface near the T point, eq.(14l. Hence, the effects obtained within this approach are 



clearly associated to the proximity to the edge of the Brillouin zone. 

In Fig. [5] we present the contribution to the LDOS calculated by allowing the tip to describe a circular path around 
a Nb atom, as a function of the angle 9$. There are clear position dependent changes, produced by the degeneracy of 
the Nb orbitals, combined with the modulation of the superconducting gap. 

The changes in the electronic structure, as measured by the STM tip, near a vortex can be obtained by solving 
the Eliashberg equations in a vector potential. They require a knowledge of the modulation of the gap, which will 
be reflected in the density of states. Qualitatively, and using the simpler Ginzburg-Landau equations, it can be 
expected that the non trivial hybridization around the Fermi surface will lead to changes in the gradient term in the 
Ginzburg-Landau expansion. This term modifies the superconducting gap when the superconducting order parameter 
is modulated spatially over distances comparable to the coherence length, as it happens near a vortex. The gap is 
reduced. A proportional reduction of A and Aj near the core of a vortex leads to changes in the density of states in 
qualitative agreement with the experiments. 



CONCLUSIONS 



In summary, we have observed a sizable modulation of the superconducting density of states at interatomic distances, 
localized at energies corresponding to the size of the different superconducting gaps found over the Fermi surface of 
2H-NbSe2, and at reciprocal space vectors fully coinciding with the atomic Se surface lattice. Moreover, the star shape 
vortex structure found by Hess et alP^l, also shows a significant atomic size modulation, not detected previously, close 
to the vortex core. These modulations appear again localized at the lattice reciprocal space, with a small contribution 
at spacings corresponding to the CDW, and at energies that are however closer to the Fermi level. We have performed 
a first simple analysis that shows that a combination of the degeneracies of the wavefunctions near high symmetry 
points, and a finite angular dependence of the superconducting gap around the Fermi surface, are enough to induce an 
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angular modulation on the current measured by an STM tip, only at voltages within the superconducting gap. Other 
approaches involving different symmetries for the tip substrate coupling, or different superconducting gap modulations 
can also lead to a spatially varying local superconducting density of states. A more sophisticated microscopic model 
taking into account more precisely the band structure of 2H-NbSe 2 and possible tip substrate coupling functions should 
lead, through the comparison to the experiment, to new, additional insight into the properties of the superconducting 
phase of 2H-NbSc2. From the present approach, it already becomes clear that the atomic size features in the STS 
curves open a new route to measure the reciprocal space structure of the superconducting gap, not only in 2H-NbSe 2 , 
but also in other systems with a significant gap anisotropy. In the present work we have identified the two energy 
scales (Vi = 0.75 mV and V 2 = 1.2 mV) corresponding to the two main gap magnitudes open over the Fermi surface 
in this compound, as well as a six fold modulation associated with them. We also give first evidence for an additional 
atomic size modulation of the superconducting gap at the vortex cores of 2H-NbSe 2 . 

The present experiment breaks with the conventionally accepted view that the intrinsic superconducting properties 
are homogeneous at length scales below the superconducting coherence length. The atomic scale structures found at 
zero field and within vortices are a surprising feature of superconductivity, which is possibly unique among macroscopic 
quantum ordered states. 
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